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We study the dynamics of a one-component liquid constrained on a spherical substrate, a 2-sphere, 
and investigate how the mode-coupling theory (MCT) can describe the new features brought by the 
presence of curvature. To this end we have derived the MCT equations in a spherical geometry. We 
find that, as seen from the MCT, the slow dynamics of liquids in curved space at low temperature 
does not qualitatively differ from that of glass-forming liquids in Euclidean space. The MCT predicts 
the right trend for the evolution of the relaxation slowdown with curvature but is dramatically off 
at a quantitative level. 


I. INTRODUCTION 

What is the influence of the curvature of space on the 
dynamics of a fluid? This question arises in several phys¬ 
ical situations, such as colloidal assemblies on droplets 
in particle-stabilized (Pickering) emulsions [1-8] or fluid 
monolayers adsorbed on large solid particles [9, 10]. Ac¬ 
cording to the frustration-based theory of glass-forming 
liquids [11, 12], depending on the dimension of space and 
the nature of the liquid, curvature may induce or release 
“geometric frustration” and therefore affect the ability to 
form glassy phases and the slowdown of relaxation. For 
spherical particles in 2 dimensions, hexagonal or hex- 
atic 6-fold order is prevalent in flat (Euclidean) space, 
and a constant nonzero curvature thwarts the long-range 
or quasi-long-range ordering by forcing in an irreducible 
density of topological defects [11-18]. At high density 
and/or low temperature the particle configurations then 
consist of a hexatic background interrupted by topolog¬ 
ical defects in the form of finite-length strings known as 
“grain boundary scars” [3-6, 13, 14, 16]. Correspond¬ 
ingly, there is a strong slowing down of the dynamics as 
one lowers the temperature or increases the density. One 
may wonder if this dynamical slowdown is similar or not 
to that of conventional glass-forming liquids? 

Here we address these questions through a study of the 
dynamics of a one-component atomic liquid embedded 
in a 2-sphere, ie., a surface of constant positive curva¬ 
ture, by Molecular Dynamics simulation and the mode¬ 
coupling theory (MCT). The MCT of glass formation was 
developed in the mid-80s in the context of the theory 
of simple liquids [19, 20]. In its canonical version the 
MCT predicts the dynamics of a system from the mere 
knowledge of the static pair correlation function or static 
structure factor. It is a mean-field-like theory of the glass 
transition which predicts a dynamical, ergodic to noner- 
godic, transition that is not observed in real glass-formers 
in finite dimensions where it is (at best) replaced by a 
crossover. Due to the absence of singularity in finite di¬ 
mensions, the predictions are not quantitatively accurate 
(Tc is typically overestimated by a factor of 2 in simple 
glass-forming systems). Nonetheless, the MCT appears 
to correctly capture the observed trends of the dynamics 
in physical situations where a significant change in the 


static structure factor is encountered, as, e.g., in attrac¬ 
tive colloids and other soft-matter systems [21-24]. 

The points we have therefore investigated within a gen¬ 
eralization of the MCT for a spherical geometry are the 
following: 

1) Is the predicted dynamical transition governed by a 
different type of singularity than that found for standard 
glass-forming liquids? For instance, it has been predicted 
that glass-forming liquids confined in disordered porous 
media may display within the MCT, depending on the 
model parameters, a continuous (A-type) transition dif¬ 
ferent from the discontinuous (B-type) one observed in 
bulk glass-formers [25]. Is a similar phenomenon encoun¬ 
tered for particles on a sphere due to specificities of the 
kinetics of the grain boundary scars moving in a random 
and slower evolving hexatic environment? 

2) Is the MCT able to predict the variation of the slow¬ 
ing down of relaxation with curvature? To make the com¬ 
parison more quantitative, we have carried out Molecular 
Dynamics simulations of a model one-component liquid 
on S‘^ with several values of the curvature, parametrized 
by the ratio of the radius of the embedding sphere over 
the particle diameter, and we have fitted the numerical 
data according to MCT predictions. 

We find that the MCT singularity on S‘^ is actually 
of the same type as that in Euclidean space. Thus, as 
seen from the MCT, there is nothing special about glass 
formation in curved space, which make liquids in curved 
space bona fide models for studying the generic features 
of the glass transition. As for the answer to the second 
question, we conclude that the MCT predicts the right 
trend, with the the slowdown becoming more marked as 
curvature is reduced. However, the effect is found to be 
much too weak when compared to the simulation data. 
This is another example [26] where the MCT is unable 
to properly amplify the small changes seen in the static 
pair structure to describe the observed strong differences 
in dynamics. 
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II. A SIMPLE GLASS-FORMING LIQUID ON A 
SPHERE 


A. Model and simulation 

We consider a monodisperse fluid of N atoms in¬ 
teracting through a pairwise Lennard-Jones potential 
and embedded either in a 2-sphere (the surface of a 3- 
dimensional “ball”) or in the Euclidean plane . The 
pair interaction potential is given by 


v{r) = 4e 
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where r is the geodesic distance between two atom cen¬ 
ters. The interaction is truncated at a conventional cut¬ 
off distance of 2.5cr. On a 2-sphere S‘^ {S‘^ usually refers 
to the 2-sphere of radius unity but we extend the no¬ 
tation to a sphere of any radius R), this interaction is 
also known as the “curved line of force” [9]. The units of 
mass, length, energy and time are m, a, e, and ^/mo^Je. 
On the surface of a sphere of radius R, the reduced den¬ 
sity is fixed by the number of particles and the ratio R/cr, 
which represents a dimensionless radius of curvature, as 


2N ( ( CF \\ / X 


In the Euclidean (flat) limit, R ^ oo, one then recovers 
the conventional reduced density, p {N/V)a^, where 
V is the volume (actually, the area) of the system. 

Eor studying the dynamics of the fluid on the sphere 
it is convenient to view each particle as a 3-dimensional 
rotator rigidly linked to the center of the sphere so that 
it is constrained to move at a fixed distance R of this 
center. The Hamiltonian is then 
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where uJi is the angular velocity of the rotator i, and the 
associated equations of motion read 


mwi = -r* X '^dr,v{rij ), 


( 4 ) 


with Vij the geodesic distance between rotators i and j. 

A Molecular Dynamics simulation can then be imple¬ 
mented via a “velocity Verlet algorithm”, with an update 
of the angular velocities for a time step At according to 
[27] 

u}i{t + At) =0Ji{t) - ^ [ri(t + At) X '^dr;v{rij){t + At) 

+ ri{t) X ^dr,v{rij){t)] 

( 5 ) 


and, by using Vj = Wi x r*, an update of the positions as 

Vi{t + At) = [1 + a{t, At)]ri{t) + At[u)i{t) + ^<^i] x r,, 

(6) 

where a{t, At) is determined by enforcing the constraints 
r|(t + At) = r?(t) = R^ [28, 29]. See Ref. [27] for more 
details. 

Eor illustration we have studied systems at a reduced 
density p = 0.92 for a range of temperature T and of 
ratio R/a. Eor the spherical geometry, initial configu¬ 
rations consist of particles placed randomly on S‘^ such 
that the distance between any pair of particles is always 
larger than 0.85cr. In the first stage of the simulation, 
the velocities are often rescaled in order that the mean 
kinetic energy becomes equal to the chosen temperature 
T. Eor the Euclidean geometry, we have used a rect¬ 
angular cell of aspect ratio equal to 2/^/3 with periodic 
boundary conditions to minimize the number of topolog¬ 
ical defects. The initial configuration is then taken as a 
triangular lattice. As for 5'^, the velocities are rescaled 
in the early stage of the simulation in order to enforce 
that the system stays at a given temperature T. 


B. Structure and dynamics 

We have characterized the structure and the dynamics 
of the system by computing two quantities, the static 
structure factor S{k) and the self-intermediate scattering 
function Fs{kR). (Other observables are considered and 
described in a previous paper [27].) In 5'^, the static 
structure factor is expressed as (see also the Appendix A 
for details) 


where k is an integer, is the kth Legendre polynomial, 
and Vij is the geodesic distance between atoms i and j. 
The self-intermediate scattering function is given by [18] 



where rj(0, t) is the geodesic distance traveled by atom j 
between times 0 and t. 

The static structure factor is the necessary input for 
the MCT equations. We have obtained it from Molecu¬ 
lar Dynamics simulation for different system sizes, N = 
1000, 2000, 4000, and 12000, on S‘^. Eor the fixed cho¬ 
sen reduced density p = 0.92 these sizes correspond to 
an increasing reduced radius of curvature, Rja 9.3, 
13.2, 18.5 and 32.3. Note that the chosen system sizes, 
with N > 1000 are large enough to avoid trivial finite- 
size effects that are unrelated to curvature and would be 
already present in the Euclidean plane. In addition we 
have also studied the one-component liquid on the Eu¬ 
clidean plane, above and near the ordering transition. In 
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FIG. 1. Static structure factors S{k) of the one-component 
Lennard-Jones liquid in 2-dimensional space versus the re¬ 
duced wavenumber kajR. (a) Effect of temperature at a con¬ 
stant curvature R/a = 18.5: T = 1, 1.3 and 1.5. (b) Effect 
of curvature at constant temperature T = 1.45: R/a = 9.3, 
R/a — 32.3 and Euclidean case (in this case the data is plot¬ 
ted versus ka where k is the wavevector modulus). In all cases 
the reduced density is p — 0.92. The insets show a zoom in 
on the first peak. 


this case we have used systems that are sufficiently large, 
up to 14400, so that the static pair correlation function 
reaches its asymptotic value for a distance less than half 
the linear size of the simulation cell at the temperature 
studied. 

In Fig. 1, we display the effect of the temperature 
at constant curvature (a) and the effect of curvature at 
constant temperature (b) on the static structure factor. 
As can be seen, the variation is modest and unremarkable 
in all cases. The main variation appears through the 
value of the first peak, especially when varying curvature 
at fixed temperature. Correspondingly, a change in the 
(split) second peak is also observed. 

We illustrate in Fig. 2 the evolution with temperature 
of the self-intermediate scattering function Fs(/c,t), for 
k corresponding to the typical interatomic distance, in 
the case R/a = 18.5. With decreasing temperature, one 



FIG. 2. Self-intermediate scattering function Fs{kR) versus 
time on a log-linear plot for R/a = 18.5. The temperature 
varies from T = 3 (left curve) down to T = 0.6 (right curve). 




FIG. 3. (a) Arrhenius plot of the relaxation time r versus 1/T 
for R/a = 18.5. (b) Variation with curvature of the relaxation 
time for a fixed temperature T = 0.9. Note the difference in 
scale for the y-axis in the two plots. 


can see the gradual change from a simple exponential 
decay to a more stretched, two-step-like, decay. How¬ 
ever, in the range of temperature under study, there is 
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no well-defined plateau in the scattering function. Fig. 
3 shows the relaxation time r as a function of 1/T for a 
given curvature, R/a = 18.5, and as a function of cur¬ 
vature for a given temperature, T = 0.9 (r is defined 
from Fs{k,T) = 0.1). As expected, curvature thwarts 
crystallization and the systems remains liquid down to 
temperatures below the ordering transition in the Eu¬ 
clidean plane (which is around 1.3 — 1.4). One also finds, 
as already observed for a liquid in a surface of constant 
negative curvature [15-17], that the relaxation time in¬ 
creases with decreasing curvature, i.e. with increasing 
R/a. 


III. MODE COUPLING THEORY IN 
SPHERICAL GEOMETRY 

A. Preliminaries 

To derive the generalization of the MCT equations for 
a liquid embedded in a 2-sphere, we repeat the steps of 
the derivation in Euclidean space and introduce the ap¬ 
propriate changes associated with the positive curvature. 

Eor a monodisperse system on 5'^, the microscopic den¬ 
sity p(fl) is defined as 

^ 1 

p{n = {0,cp)) = j2^m-^i), 

i=l 

= (9) 


47r/A (Yfc o(^i)Ffc o(^j)) [which coincides with Eq. 

( 7 )]. 

The MCT starts with a generalized Langevin equa¬ 
tion for the Eourier components of the density and 
their time derivative. One first defines the vector 
^Pk,i 

trix = (Al Because of the space 

homogeneity, Ck,i{t) = Ck,o{t) = Ck{t), and 


j and the associated correlation ma- 



Ck{t) = 


NF{k,t) 
2/c-|-1 

N ■ 


N 


{SploiO)Sh,o{t)) 


(13) 


where a dot indicates a time derivative. 

By using the projection-operator technique, the evolu¬ 
tion equation for the correlation matrix can be cast in the 
form of a generalized Langevin equation [30], from which 
it is easy to derive the following equation for the normal¬ 
ized dynamic structure factor /(/c,t) = F(kA)/S{k)\ 




AnksT k{k + 1) 
mV S{k) 


f{k,t)+ 



,t- t')f{k,t') = 0 


(14) 


where V = is the surface area of the sphere, the 

memory function K(kA) is defined as 




and the so-called “random force” R/c(t) is given by 


where 0 and (j) are the azimuthal and polar angles, respec¬ 
tively. By using Eq. (Al), the Eourier transform pk^i of 
the density reads 


I —z— ^ 

= (10) 

i—\ 

where k and I are integers such that k > 0 and —k< 
I < k, and Yj/i is (the complex conjugate of) a spherical 
harmonics (see Appendix A). The dynamic structure 
factor is defined from the density fluctuations Sp = p—{p) 
as 


= V E (11) 

^ l=-k 

where the bracket denotes the thermal average. Because 
of the homogeneity of the system, (^6ph{0)6pk,i{t)) = 

{^Pk,oiO)^Pk,o{t)) and 

F{k,t) = ^^{Splo{0)Spk,o{t))- ( 12 ) 

At time t = 0 the above function reduces to the 
static structure factor, F(kA = 0) = S{k) = 


-Rfc(f) (16) 


with V the projection operator onto A/c(0), iC the Liou- 
ville operator, and 


Rk{^) 


\/2fc + 1 


Spk,o{0) 


4TTkBTk{k + l) 
mV S{k) 


(5pfe,o(0) 


(17) 


B. Mode-coupling approximation 

The above generalized Langevin equation is exact but 
purely formal. The MCT of glass-forming liquids consists 
of the following approximations: (i) One only retains in 
the dynamics of the random force the contribution of 
bilinear combinations of the density fluctuations; (ii) the 
operator is replaced by V 2 F^^V 2 where V 2 is a 

projection operator defined for any variable X as 

E E Sk.k. (B£.k.X) 

ki7^k2* k37^k4* 

(18) 

where the vector B has components ^kika = ^Pki^Pka 
(we have simplified the notations by writing k = (/c,/), 
where /c > 0, |/| < /c, and its “conjugate” k* = (/c, — /). ). 
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(iii) Finally, the four-point density correlation function is 
factorized into a product of two-point correlation func¬ 
tions. (Note that in principle the output also involves the 
static triplet correlation function but the latter is almost 
always considered in a factorized approximation.) 

After some lengthy but straightforward algebraic ma¬ 
nipulations that are detailed in Appendix B, we have 
obtained the following expression for the mode-coupling 
memory function on S‘^: 


K{k,t) 


nphsT 

2mV‘^k{k + 1) 


k-\-k\ 

E E Zk{kiM)f{ki,t)f{k2,t) 

fel>0/c2 = |fe-fel| 


(19) 


with 


Zk{k^,k2) = (2fci + l){2k2 + 1) X 

(^k{k + 1) + ki{ki -|- 1) — ^ 2(^2 + l))c(A)i)+ 


(^k{k + 1) + k2(k2 + 1) 


^ 1(^1 + 1))c(A:2) 


S{k^)S{k2). 

( 20 ) 


In the above expression denotes a Wigner 

3j symbol and c{k) is the direct correlation function 
which is related to the static structure factor S{k) by 
the Ornstein-Zernike equation 1/S{k) = 1 — pc{k). Note 
that the Ornstein-Zernike equation takes the same form 
as in the Euclidean case thanks to the appropriate choice 
of the Fourier transform in S‘^. 

An equation can also be derived for the self¬ 
intermediate scattering function Fs{k,t) [20], which is 
more conveniently computed in simulations. However, 
it is known that for wavenumbers corresponding to typi¬ 
cal interatomic distances the temperature dependences of 
the relaxation times associated with the dynamic struc¬ 
ture factor and the self-intermediate scattering function 
are similar. More importantly for our purpose here, the 
critical temperatures Tc at which the dynamics are pre¬ 
dicted to freeze are identical. 


C. Euclidean limit 

In the two-dimensional euclidean space the MCT equa¬ 
tions have already been studied [31, 32] and read 

f{k,t) + ^^^-^rj-f{k,t) + j dt'K{k,t-t')f{k,t') =0 
m D[K) Jo 

( 21 ) 

with 

f{ki,t)f{\ki-k\,t) (22) 


and 


xV k.ki , k.(k — ki) 

^ki,k-ki ^ k 

To facilitate the comparison of the above expressions with 
the limit of the spherical case derived in the previous 
subsection when ^ 00 , we rewrite the function K{k^ t) 
as 


K{k,t) = /y dk,k, dk2k2Vl^,^ 

X Sik^)sik2)fik^,t)fik2,t) ^24) 

^y-(k^+kf+k^)+2(k‘^kl+k‘^kl+klkl) ^ ^ 

with a new vertex Vk^^k 2 defined as 

Vk,M = (k^ -kl + kl)c{k^) + (fc2 + kl - kl)c{k2). (25) 

Note that Eqs. (24) and (25) only depend on the modulus 
of the wavevectors. 

When R ^ 00 in Eqs. (14), (19) and (20), the 
wavenumbers /c, /ci, • • •, which all go to infinity, must be 
rescaled by i?, be., k = Rk with k G Under such 
a rescaling, the structure factor on 5'^, S{k) converges 
to S{k), where the latter is now the Euclidean structure 
factor, and similarly F{k, t) converges to its Euclidean 
counterpart F{k, t). The coefficient of the term in /(/c, t) 
in Eq. (14), which is like the square of an effective fre¬ 
quency, goes to ksTk^/[mS{k)], which is precisely the 
coefficient in the Euclidean equation, Eq. (21). 

To treat the memory term, one must use the large- 
wavenumber limit of the Wigner 3j symbols, [33] 

Al k2 ksV 
\0 0 0 ) 

~- , ^ (26) 

T^V-ikf + q + kj) + 2{klkl + klkl + k^k^) 

where \ki — k 2 \ < ks < ki F k 2 and ki, k 2 ^ ks ^ 00 . The 
3j symbols in Eq. (26) are only nonzero when the sum 
kl F k 2 F ks is even. 

After introducing the above property in Eq. (20), using 
the rescaling of the wavenumbers, and replacing the sums 
in Eq. (19) by integrals, be.. 


E E 

kl >0 k2 = \^~^i I 


k-\-ki / 7 N ^ 

kl k2 ks 


0 0 0 
2 


-I nOO nk-\-ki 

2 / /- - 

^ Jo J\k — ki I 


dk2 


T^\j-{ki + H + F + Flkl + kF + kF) 


■ ■ , (27) 


one arrives at Eqs. (24,25). Note that there is a fac¬ 
tor 1 /2 in front of the integrals because of the condition 
that kl F k 2 F ks is even. This completes the proof that 
the MCT equations on 5'^, which we have derived above, 
converge to the MCT equations of the Euclidean plane 
when the radius R of the sphere goes to infinity. 
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IV. MODE-COUPLING SINGULARITY AND 
COMPARISON WITH SIMULATION DATA 


A. Predicted dynamical singularity 


As mentioned in the Introduction, our goal is two¬ 
fold. On the one hand, we would like to study whether 
curvature induces a slowdown of relaxation whose na¬ 
ture is qualitatively different from that of standard glass¬ 
forming liquids. Such a qualitative difference could be 
signaled by a dynamical mode-coupling transition of a 
different character than the discontinuous one predicted 
for glass-forming liquids in 2- and 3-dimensional Eu¬ 
clidean spaces. On the other hand, increasing curvature 
is found to speed up the dynamics and we want to check 
if, and to which extent, the MCT is able to capture this 
trend. For studying these questions, it is sufficient to 
consider the nature and the curvature dependence of the 
dynamical singularity of the MCT equations. 

We therefore consider the long-time limit f{k) = 
limt^oo of the MCT equation, which is known as 

the nonergodicity parameter. Within the MCT frame¬ 
work [20? ], a dynamical transition at some critical Tc 
separates a high-temperature, ergodic phase, in which 
f{k) is equal to zero, from a low-temperature, noner- 
godic phase in which f{k) > 0. For glass-forming liquids 
in Euclidean space, the transition is discontinuous, with 
a jump of f{k) at Tc (a B-type singularity), whereas in 
other situations, such as liquids confined in disordered 
environments, the transition may be continuous (an A- 
type singularity) [25]. 

The equation for the nonergodicity parameter is easily 
obtained from the full MCT equation in the long-time 
limit by taking the time Laplace transform, and one ob¬ 
tains that [20] 

(28) 

where, for a liquid on 5'^, 

= SVkHkh)^ F Zk{kuk2)m)f{k2) (29) 

^ ^ ki,k2 


with Zk{ki^k 2 ) given in Eq. (20). 

A numerical solution of Eqs. (28,29) can be obtained 
with an iterative procedure. 


1 -/(*+!) (fc) 


Mg\k), 


where M^{k) is calculated by using the nonergodic¬ 
ity parameter of the previous iteration {f^^\k')}. For 
the sake of simplicity, the initial condition is chosen as 
f^^\k) = 1. The iterative procedure is stopped when 
the condition — f^'^\k)\‘^ < 10“^ is fulfilled. 

The convergence is usually obtained after few iterations 
except close to the transition where the number of itera¬ 
tions is of order 100. 



T 


FIG. 4. MCT prediction for the temperature dependence of 
the nonergodicity parameter /max = /(A^max) for S‘^ with dif¬ 
ferent values of R/a and for the Euclidean case. For the 
largest sphere, R/a = 32.3, the numerical uncertainty is 
larger which explains that the magnitude of the discontinuity 
at Tc is slightly too small. 


For the Euclidean model, the memory term is given by 

Moo(fc) = ^ r flt-th dk2k2Vl,^ 

S(ki)S(k2)f(ki)f(k2) /3Q^ 

■y/-(k^ + kf + k^)+2(k^kf+k2k^+kfk^) ‘ ^ ^ 

The numerical integration is a little more subtle, be¬ 
cause at a given ki the integrand has singularities at the 
bounds of the integral over k 2 . By using a trapezoidal 
rule, the numerical method goes to the exact result as 
the inverse of the square root of the number of points, 
which cannot provide a sufficient accuracy for the itera¬ 
tive method. To improve the computation of the integral, 
we have estimated it by using two different steps for the 
integration, and it is then easy to show that one can ob¬ 
tain a convergence to the exact value as the number of 
chosen points to the power 3/2. 

The solution of Eqs. (28), (29), and (30) is shown in 
Fig. 4 where we plot the temperature dependence of the 
nonergodicity parameter /^ax = /(^max), where fc^ax 
corresponds to the value of the first peak of S{k), for 
several values of the reduced radius of curvature and for 
the Euclidean case. Note that in the latter case, the pre¬ 
dicted dynamical transition takes place at a temperature 
that is higher than the ordering transition in the plane 
so that the associated static structure factors are those 
of a liquid. (We have varied the system size and checked 
that the finite-size effects are negligible for the systems 
under study.) 

The first observation is that the transition is discon¬ 
tinuous in S‘^ just as it is in Euclidean space. Therefore, 
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from the MCT perspective, there is nothing anomalous 
in the glass transition of liquids in curved space. The sec¬ 
ond observation is that the predicted critical temperature 
Tc{R) decreases with increasing curvature, from 1.49 in 
the Euclidean plane to 1.24 for the largest curvature cor¬ 
responding to R/a = 9.3. In Fig. 5 we also illustrate 
the /^-dependence of the nonergodicity parameter at the 
transition Tc for the Euclidean plane and one curvature. 



FIG. 5. MCT prediction of the nonergodicity parameter f(k) 
versus the reduced wave number (ka/R for and ka for E^) 
at the transition T = T“, for with R/a = 9.3 and 32.3, 
and for the Euclidean case. 

From the above MCT equations, one can derive some 
analytical results. First, as the dynamical singularity in 
is of the same B-type as the singularity found for 
glass-forming liquids in Euclidean space, one finds the 
same critical behavior with in particular a square root 
dependence of the nonergodicity parameter in the glass 
phase, 

(31) 

when T T ~. 

In addition, one can study the asymptotic behavior of 
the MCT equations when R ^ oc. We find that the first 
correction to the kernel K{k, t) is in cr/R^ which leads to 
an expected behavior of the transition temperature as 

T,(oo)-T,(i?)~^ (32) 

when approaching the flat (Euclidean) space limit. We 
plot in Fig. 6 Tc(oc) —Tc{R) versus a/R. Within the nu¬ 
merical accuracy the results are compatible with a linear 
behavior. 


B. Comparison with simulation data 

We would like to compare the predictions of the MCT 
with the simulation data. This is not an easy task as 



FIG. 6. Gurvature dependence of the predicted MGT critical 
temperature: Tc{oo) —Tc{R) versus a/R for R/a = 9.3, 13.2, 
18.5 and 32.3. The value of the transition temperature on the 
Euclidean plane is Tc (oc) ~ 1.49. The full line is the expected 
linear asymptotic behavior. 


the dynamical transition predicted by the MCT is not 
observed in finite-dimensional systems and, as a result, 
the predictions are quantitatively inaccurate. The con¬ 
ventional procedure to compare theory and simulation is 
to try to empirically fit the numerical data to the the¬ 
oretical results by letting the temperature of the postu¬ 
lated dynamical singularity be an adjustable parameter. 
This empirical temperature, , can be further used 

as a rescaling parameter to improve the quality of the 
comparison. Our interest here is not in pursuing such a 
rescaling, but the empirically determined is a use¬ 

ful quantity to compare predicted and observed trends in 
a restricted domain of relaxation times. 

From the simulation data for the intermediate scat¬ 
tering function we have extracted an estimate of the re¬ 
laxation time r. We have then fitted the temperature 
dependence of the latter to the MCT prediction by us¬ 
ing the scaling law r ex (T — Tc)~^ over an “optimal” 
range of relaxation time. To restrict the number of ad¬ 
justable parameters, we have set 7 = 2.38 according to 
results obtained in the two-dimensional Euclidean space 
[31]. Tc = is then obtained from the best fit to 

the data. 

Fig. 7a displays a log-log plot of the simulation data 
with the associated MCT fits for the 4 different curva¬ 
tures studied. (One observes that the prefactor of the 
scaling law depends on the reduced curvature and de¬ 
creases with R/(7 but the physical significance of this 
trend is unclear.) In addition. Fig. 7b illustrates the 
same MCT fit now shown on an Arrhenius plot for one 
curvature. As well known, the quality of the fit is rather 
poor and is only valid over at most a couple of orders of 
magnitude in r [26]. Nothing here is specific of curved 
space. 

The predicted and empirically determined critical tern- 











(T-T^)/T^ 



FIG. 7. a) Log-log plot of the relaxation time r versus (T — 
rjn(emp)yrjniemp) different curvatures. is chosen to 

obtain the best collapse of the simulation data to the scaling 
law T (X {T — Tc)~^ with the value of the exponent hxed to 
7 = 2.38. b) Illustration of the same MCT ht (full line) on 
an Arrhenius plot of log(T) versus 1/T for R/a = 18.5. 


but the quantitative estimate of the effect is completely 
off. 


V. CONCLUSION 

We have studied the dynamics of a one-component liq¬ 
uid constrained on a spherical substrate and investigated 
how the mode-coupling theory (MCT) can describe the 
new features brought by the presence of curvature. To 
this end we have derived the MCT equations for the time- 
dependent pair correlation function of the density fluctu¬ 
ations in a spherical geometry. The ratio of the radius 
of the embedding 2-sphere over the atomic diameter is a 
new control parameter. 

We find that, as seen from the MCT, the slow dynam¬ 
ics of liquids in curved space at low temperature does not 
qualitatively differ from that of glass-forming liquids in 
Euclidean space [34]. The dynamical transition is of the 
discontinuous (B) type in all cases. The MCT predicts 
the correct trend for the evolution of the relaxation slow¬ 
down with curvature but is dramatically off at a quanti¬ 
tative level. As found in other cases [26], the MCT is un¬ 
able to describe situations where significant differences in 
the temperature dependence of the dynamics come with 
only modest changes in the static pair structure. 
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Appendix A: Fourier transform on 

The sphere is a manifold with a finite area, so that 
the Fourier transform of a function of the angular coor¬ 
dinates, /(O), is given by a discrete number 


peratures are given in Table 1 for the 4 different curva¬ 
tures under study. Several comments can be made: 

1) As already found in all other MCT studies of glass¬ 
forming liquids, the theoretically determined Tc is much 
higher than the empirically determined one, with a factor 
2 or more between them. 

2) The theory predicts the correct trend with increas¬ 
ing curvature, z.e., that the slowdown becomes weaker 
as curvature increases: this is signaled by a decrease of 
Tc{R) with decreasing R/a. 

3) The magnitude of the variation with curvature is 
dramatically underestimated by the theory: the reduc¬ 
tion in Tc between the lowest nonzero curvature {R/a = 
32.3) and the highest one {R/a = 9.3) is about 4 times 
too small in the theory. The qualitative trend is right 


fk,i dny7(n)/(n), (Ai) 

where k and I are integers such that k > 0 and —k< 
I < k, and Y^i{Q) denotes the complex conjugate of the 


TABLE 1. Table of temperatures 


rji rp{emp) rji jrjn{emp) 


9.3 

1.24 

0.45 

2.8 

13.2 

1.33 

0.53 

2.7 

18.5 

1.38 

0.62 

2.5 

32.3 

1.43 

0.70 

2.2 
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spherical harmonics Yk^i{Q) which is defined as 


The MCT memory function is obtained from 


(A2) 

where Pk^i{cosO) is an associated Legendre function of 
first kind. With this definition, lo,o = 1 /V^ and 

f^2 

The corresponding inverse Fourier transform is then 

I - k 

/(**) =vE E Yu(a)hj. (A 3 ) 

k>0l=-k 

For an isotropic function f{Q) = /(6>), fk,i = fk,o = 
f{k). In the limit ^ oo, after having introduced k = 
k/R and /(r) = f{r/R), one obtains that f{k) goes to 

POO 

f{k) = 2 'k drrJo{kr)f{r), (A 4 ) 

Jo 

which is the usual Fourier transform in 2 dimensions. A 
similar reasoning holds for recovering the inverse Fourier 
transform. 


Appendix B: Mode-coupling approximation on 

We summarize here the main steps of the derivation of 
the MCT equations on S‘^. 

The projection operator V2 on the pair products of 
density modes, when applied to the random force, leads 
to 


r 2 Rk= U(ki,k2)Bk,k,, (Bl) 

ki7^k2* 


with 


U(ki,k2)= ^ 

k3 7^k4* 

(B 2 ) 

where [-D]kik2,k3k4 = (-Bkik2-®k3k4) and we recall the 
notation = {ki, k) and k* = {ki, —k)- 

Within the Gaussian approximation, the 4 -point den¬ 
sity correlation function [T^]kik2,k3k4 5 with ki 7^ k2* and 
ks / k4*, is given by 


[L^]kik2,k3k4 ~ (^Pki^Pk3)(^Pk2^/^k4) 

+{SptJp^,){SplJpu,) 


— ('^ki,k3^k2,k4 “t“'^ki,k4^k2,k3) 

2 S{k,)S{k2) 

{2ki + l){2k2 + 1) ’ 


(B 3 ) 


so that 


[D- 


Jkik2,k3k4 — 


{2ki + 1)(2A:2 + 1) 
4N^S{ki)S{k2) 


('^ki,k3^k2,k4 "1“ '^ki,k4^k2,k3 


(B 4 ) 


which, by using the above expressions and, once more, 
the Gaussian approximation, can be rewritten as 

(B6) 

Given the definition of the random force in Eq. ( 17 ), 
the calculation of (^Pk3^Pk4^^) appearing in Eq. (B 2 ) 
requires the consideration of two terms: 

C/k 3 ,k 4 ,k = (SpUp^Jp^), (B 7 ) 

rk3.k4,k = (B 8 ) 

where, here and below, k = (/c, 0). 

The 3 -point density correlation function Tk3,k4,k is ex¬ 
pressed in the convolution approximation as 

(B 9 ) 

The other correlator I/k3,k4,k can be rewritten by using 
the definition of the Eourier transform for the density 
(see Appendix A) and the following expression for its 
derivative: 



(BIO) 


^k3,k4,k is then expressed in terms of sums over i and 
j of (VTk3(^i)VTk(^i)hk4(%)). This average is easily 
calculated thanks to the three following properties: 

1) A product of gradients is expressed with only Lapla- 
cians through 

Vf.Vg=^[Aifg)-fAg-gAf]. (Bll) 

2) The spherical harmonics are Eigenfunctions of the 
Laplacian: AYk^i = —k{k + 

3 ) A product of two spherical harmonics is a linear 
combination of spherical harmonics with coefficients de¬ 
pending on 3 j symbols: 


Yk,iYk^^i^ = 



(2fc + l)( 2 fc 3 + l)(2fc' + l) 
Air 

ks k k\ 

0 0 0 / 


(B 12 ) 


One then obtains 

t^k3,k4,k = ^^^ 2TT[k{k + 1) + k3{k3 + 1) - k^iki + 1)] 

/ \ 

ks k 4 ^ k 



0 0 0 


S{k 4 ). (B 13 ) 
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Summing up all the preceding results and introducing After inserting the above result in Eq. (B6), one readily 

the direct correlation function c{k) = (l/p)[l — 1/S{k)] finds the expressions in Eqs. (19,20) of the main text, 

then lead to 


V/c(ki,k2) — ^ Tr V(2/c + l)(2/ci + 1)(2 /c2 + 1)X 

mV 

^2(^2 “1“ l)]c(A:i) + [k{k + 1) + k 2 {k 2 + 1) — ki{ki + l)]c(A:2)^. 

(B14) 
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